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LONGITUDINAL  RESISTIVITY  AND  INTERMEDIATE  MACH  NUMBER 
COLLISIONLESS  TRANSVERSE  MAGNETO  SO  NIC  SHOCKS 

I.  Introduction 

It  is  now  reasonably  well  established  that  transverse  collisionless 

shocks  can  exist  in  at  least  two  Mach  number  regimes.  For  m  (the 

magnetosonic  Mach  number)  less  than  about  2.5,  transverse  resistivity  is 

the  dissipation  mechanism  which  forms  the  shock. Since  the  plasma  is 

nearly  collisionless,  this  resistivity  is  anomalous  in  that  it  is  assumed 

to  arise  from  an  electron-ion  streaming  instability  in  the  shock  front. 

Measurements  of  the  turbulent  spectrum  in  the  shock  front  show  that 

fluctuations  with  wavenumber  k  ~  kp,  the  Debye  wavenumber,  are  indeed 

excited.  ‘  It  is  generally  assumed  that  the  resistive  dissipation  heats 

only  the  electrons.  However  a  small  amount  of  ion  heating  cannot  be  ruled 

out.  It  is  certainly  not  easy  to  see  how  an  instability  can  heat  only  one 

species.  An  ion  acoustic  instability,  at  least,  does  give  rise  to  a  small 

amount  of  ion  heating.^  In  fact  one  theory  of  the  resistive  shock 
o 

structure  has  this  small  ion  heating  playing  an  important  role  in 
determining  shock  structure. 

For  larger  m,  it  is  now  well  established  both  from  theory®’ *®* 11  and 
12—19 

particle  simulations,  4  that  ion  reflection  from  the  shock  front  can  act 

as  a  dissipation  mechanism  in  one  of  two  ways.  On  the  one  hand,  these 

reflected  ions  may  stream  ahead  of  the  shock  if  they  are  unmagnetized  or 

sufficiently  weakly  magnetized.®* ^^“15  if  this  be  the  case,  the 

standard  Rankine-Hugoniot  conservation  relations  do  not  apply  because  of 

the  additional  species  upstream.  On  the  other  hand,  the  reflected  ions  may 

turn  around  in  the  magnetic  field  and  drift  downstream,  their  large  Larmor 

orbit  velocity  effectively  becoming  part  of  the  downstream 

pressure. »16-19  Since  there  is  no  additional  upstream  species  in  this 
Manuscript  approved  October  27,  1983. 
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case,  Rankine-Hugoniot  conservation  relations  do  apply.  In  fact  the 
simulations  of  Refs.  18-19  were  successful  because  the  Rankine-Hugoniot 
relations  were  exploited  in  initializing  the  system.  Of  course  for  the 
latter  dissipation  mechanism  to  be  significant,  the  time  of  the  experiment 
or  simulation  must  be  longer  than  about  ir/u)cj,  where  is  the  ion 
cyclotron  frequency. 

The  purpose  of  this  paper  is  to  show  that  between  these  two  ranges  of 

Mach  number,  there  exists  an  intermediate  range.  In  this  intermediate 

range,  the  shock  structure  and  dissipation  are  dominated  by  the  breakdown 

of  the  quasi-neutral  approximation.  Above  a  critical  Mach  number,  the 

quasi-neutral  fluid  equations  predict  infinite  slope,  implying,  of  course, 

the  breakdown  of  quasi-neutrality.  However,  the  full  set  of  fluid 

equations,  including  charge  separation,  are  nonsingular.  Thus,  above  this 

Mach  number,  the  flow  sets  up  an  electrostatic  ion  acoustic  wave  train. 

The  Landau  damping  of  this  wave  is  the  dissipation  mechanism  which  forms 

the  shock.  This  Landau  damping  may  heat  either  the  electrons  or  ions.  In 

purely  electrostatic  shocks  with  very  cold  ions,  the  electron  Landau 

20 

damping  seems  to  dominate.  Since  the  electron  Landau  damping  rate  is 

small,  the  shock  structure  is  oscilliatory.  On  the  other  hand,  in 

magnetosonic  shocks  the  ions  are  warmer,  so  ion  Landau  damping  dominates 
and  the  shock  structure  is  monotonic.  It  is  interesting  that  this 
dissipation  mechanism  does  not  rely  on  exciting  some  plasma  instability  in 
the  shock  front,  and  in  that  sense  it  is  not  anomalous.  Indeed,  the 
electrostatic  subshock  has  a  Debye  length  scale  size,  so  it  is  almost 
impossible  that  an  instability  with  still  smaller  wavelength  could  be 
excited. 


While  this  dissipation  mechanism  is  not  anomalous  it  is  very 
complicated*  The  simplest  model  of  it  shows  tnat  it  is  nonlocal,  and  in 
that  sense  non  fluid  like.  However,  there  is  a  new  simple  fluid 
dissipation  model.  When  the  electrostatic  oscillation  is  set  up,  the 
longitudinal  electron  and  ion  velocities  are  no  longer  equal.  A 
longitudinal  friction  force  between  the  two  is  an  additional  possible 
dissipation  mechanism  which  apparently  has  not  been  examined  in  the 
previous  literature  on  high  Mach  shock  waves.  An  important  purpose  of  this 
paper  is  to  investigate  the  effect  of  such  longitudinal  resistivity  on 
shock  structures.  Our  conclusion  is  that  it  is  a  very  good  fluid  model  for 
collisionless  magnetosonic  shocks  with  Mach  number  in  the  intermediate 
regime.  It  has  a  precise  onset,  namely  when  the  shock  reaches  the  critical 
Mach  number.  It  shows  that  shocks  above  this  Mach  number  have  an 
electrostatic  subshock  embedded  within  in  and  its  structure  is  governed  by 
longitudinal  resistivity.  When  the  Mach  number  increases  to  a  second 
critical  value,  the  full  fluid  equations  become  singular,  and  this  defines 
the  upper  limit  of  the  intermediate  range.  The  interpretation  of  this 
second  singularity  is  simply  made  in  terms  of  ion  reflection.  Above  this 
point,  the  shock,  if  it  exists,  is  governed  by  ion  reflection  as  described 
in,  for  Instance,  Refs.  9-19.  At  no  point  is  it  either  necessary  or 
physically  correct  to  postulate  a  viscosity. 

While  the  theory  we  deal  with  here  concerns  only  steady  state  shocks, 
the  general  conclusions  should  also  apply  to  non-steady  flow  of 
Interpenetrating  plasmas.  Specifically,  electrostatic  interactions  should 

be  Important  in  a  recent  series  of  laboratory  experiments  on  the  flow  of  a 

2 1 

laser  produced  plasma  through  a  background  magnetized  plasma. 
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In  Section  II  we  review  experimental  data  on  high  Mach  shock 

structures  and  show  that  it  indicates  the  existence  of  this  intermediate 

regime*  To  do  so,  we  note  that  in  this  Mach  number  range  many  of  the 

experiments  are  over  in  a  time  of  order  of  or  less  than  ti/oj  .  ,  so 

cl 

1  3  22  23 

convection  of  reflected  ions  downstream  *  »»  cannot  be  the  dissipation 
mechanism  of  importance*  Also,  downstream  ion  temperatures  are  measured  in 

several  experiments  and  they  agree  with  Rankine-Hugoniot.  This  means  that 

3  22 

pure  reflection  of  unmagnetized  ions,  '  cannot  be  the  dissipation 

mechanism  either. 

In  Section  III  we  discuss  shock  structure  for  unmagnetized  ions  where 
the  quasi-neutrality  assumption  is  valid  and  also  where  it  is  not.  We  find 
that  for  sufficiently  weakly  magnetized  plasmas,  characteristic  of  the 
experiments  in  Refs.  1-6,  22,  23,  the  electrostatic  and  magnetic  structures 
essentially  decouple,  and  that  the  electrostatic  subshock  must  ride  on  the 
crest  of  the  magnetic  shock  structure.  This  also  agrees  with  experimental 
results. 

In  Section  IV,  and  also  in  Appendix  B,  we  discuss  the  actual 
dissipation  mechanism  (Landau  damping)  and  fluid  models  for  it.  As  we 
show,  even  the  simplest  approximation  to  it  is  quite  complex,  and  does  not 
easily  fit  into  a  simple  fluid  model.  However  the  actual  physical  model 
provides  Important  insight  on  a  best  guess  for  longitudinal  resistivity  and 
on  the  partition  of  resistive  heating  between  electrons  and  ions.  Another 
interesting  phenomena  is  that  Landau  damping  can  only  be  a  valid  long  time 
description  of  the  dissipation  mechanism  if  the  shock  structure  fluctuates 
in  time.  Thus  the  shock  cannot  be  time  independent  but,  of  course,  its 
average  profile  can  be. 


Section  V  presents  the  results  of  shock  structure  calculations  in  the 
intermediate  Mach  number  regime.  Also,  where  possible,  comparisons  are 
made  with  experimental  results.  Finally,  Seccion  VI  draws  conclusions. 


In  this  Section  we  review  briefly  previous  experiments  of  transverse 
shocks  in  laboratory  plasmas.  °  >*■*»* J  Our  conclusion,  based  on 
examination  of  this  data,  is  that  for  m  <  2.5,  the  shocks  are  resistive,  as 
described  by  authors  of  these  references.  However,  there  is  a 
range  2.5  <  m  <  5  where  steady  shocks  form  which  are  above  the  critical 
Mach  number  for  resistive  shocks.  However,  the  process  of  ion  reflection 
followed  by  downstream  convection  by  the  E  x  B  drift  cannot  be  operative 
either,  if  only  because  the  duration  of  the  experiment  is  less  than  the  ion 
Larmor  time.  Furthermore,  the  dissipation  mechanism  cannot  be  simple  ion 
reflection  of  magnetized  ions  either.  For  ion  reflection  in  an 
unmagnetized  plasma,  the  second  ion  species  upstream  drastically  alters  the 
Rankine-Hugoniot  conditions  downstream.^  However,  in  two  of  the 

o  22 

experiments  cited,  downstream  ion  temperatures  were  measured  and  found 
to  agree  well  with  Rankine-Hugoniot  conditions.  Thus  our  conclusion  is 
that  there  must  be  some  sort  of  new  dissipation  mechanism  associated  with 
this  intermediate  range  of  Mach  numbers. 

In  Ref.  1,  Paul  et  al.  launch  shocks  with  Mach  numbers  2.5,  3.7  and 
6.3.  For  the  first,  resistivity  forms  the  shock  and  gives  the  electron 
heating  required  to  satisfy  the  Rankine-Hugoniot  relations.  In  the  case 
of  m  •  3.7,  the  downstream  electron  temperature  is  less  than  that  predicted 
by  Rankine-Hugoniot,  implying  that  the  ions  are  heated.  In  addition  to  the 
normal  magnetic  structure,  there  is  a  very  low  amplitude  magnetic  pedestal 
in  front  of  the  shock.  This  pedestal  is  generally  interpreted  to  mean  a 
small  number  of  ions  are  reflected  from  the  shock  front  and  convect  the 


magnetic  field  upstream.  However  in  Paul's  experiment,  the  convection  of 
these  ions  downstream  by  the  E_  x  B_  drift  is  almost  certainly  not  the 
mechanism  which  forms  the  shock  at  m  *  3.7.  The  shock  velocity  is  about 
2.5  x  107  cm/sec  and  the  shock  propagates  ahead  of  the  piston  a  distance  of 
10-15  cm,  implying  a  time  of  about  400  nsec.  However,  ir/u  ^  is  also  about 
400  n  sec.  Thus  any  reflected  ions  cannot  make  it  back  behind  the  shock 
until  the  shock  hits  the  cylindrical  axis.  However,  shocks  are  observed 
for  much  shorter  times.  For  the  measurement  at  M  *  6.3,  the  upstream 
pedestal  broadens  and  increases  in  amplitude  so  that  it  encompasses  the 
entire  shock.  At  this  Mach  number,  reflected  ions  are  almost  certainly 
playing  a  role. 

In  another  series  of  experiments,  Keilhacker  et  al.  measured  the 
properties  of  magnetosonic  shock  waves  with  Mach  numbers  between  1.5  and 
5.  Here  ion  temperatures  behind  the  shock  are  measured  by  collective 
scattering  of  laser  light  at  wavelengths  less  than  the  Debye  length.  The 
measured  ion  temperatures  agree  with  Rankine-Hugoniot  conditions.  In  this 
experiment,  the  total  transit  time  for  the  shock  is  about  ly  sec.  Thus  the 
dissipation  mechanism  cannot  be  the  convection  of  reflected  ions  downstream 
because  there  is  not  enough  time;  neither  can  it  be  the  pure  reflection  of 
unmagnetized  ions,  because  the  downstream  ion  heating  corresponds  to  that 

predicted  by  Rankine-Hugoniot. 

22 

Segre  and  Martone  have  also  investigated  transverse  shocks  in  the 
Mach  number  range  3.1  <  m  <  4.25.  In  their  experiment,  there  is  no 
magnetic  field  pedestal  leading  the  shock,  so  reflected  ions  almost 
certainly  play  no  role.  Even  if  there  were  reflected  ions,  the  ion  period 
■<*  450  n  sec,  the  propagation  time  is  only  about  250  n  sec,  so  no 
reflected  ions  could  make  it  back  downstream.  In  Ref.  22  the  downstream 
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electron  density  and  temperature  are  measured  by  laser  scattering.  However 
magnetic  measurements  are  used  to  measure  n(Te  +  T^)  downstream,  so  that 
downstream  ion  temperatures  can  be  directly  inferred.  The  measured  ion 
temperatures  do  in  fact  agree  with  that  predicted  by  Rankine-Hugoniot . 

23 

Another  series  of  shock  experiments  are  those  of  Eselevich  et  al. 

The  shocks  are  produced  in  a  tube  of  8  cm  radius  in  H,  He  or  A.  Since  the 

13  15  -3 

gas  can  be  varied,  and  the  density  varies  between  10  <  n  <  10  cm  and 

100  G  <  B  <  2000  G,  there  was  a  very  large  parameter  space  studied. 
Eselevich  et  al.  found  that  for  Mach  numbers  between  3  and  5,  steady  she  s 
form.  In  this  range  of  Mach  numbers,  a  pedestal  in  magnetic  field  might 
might  not  lead  the  shock.  For  instance,  Fig.  2a  of  Eselevich  et  al.  si 
a  Mach  4  shock  in  helium  with  no  pedestal,  whereas  Fig.  4  shows  a  Mach  4 
shock  in  hydrogen  with  a  small  (i.e.,  less  distinct  than  in  Paul's1 
experiment)  magnetic  foot  leading  the  shock.  Thus,  reflected  ions  are  not 
a  universal,  and  probably  not  a  very  important,  aspect  of  steady  transverse 
shocks  in  this  intermediate  range  of  Mach  numbers.  However,  Eselevich  et 
al. ,  like  Paul  et  al.,  did  observe  that  as  the  Mach  number  increases 
(beyond  about  5.5  generally  in  Ref.  23),  the  magnetic  foot  does  expand  out 
and  dominate  the  entire  structure. 

The  most  distinct  feature  of  Eselevich's  observations  are  that  in  the 
Intermediate  range  of  Mach  numbers,  an  isoraagnetic  potential  jump  forms  on 
the  crest  of  the  shock  wave.  The  thickness  of  the  potential  jump  is 
typically  a  Debye  length.  Thus,  a  key  feature  in  this  intermediate  range 
is  the  breakdown  of  the  assumption  of  quasi-neutrality.  The  potential  jump 
across  this  isomagnetic  transition  region  initially  increases  with  Mach 
number,  but  then  above  some  higher  Mach  number  begins  to  decrease  until 
above  a  second  critical  Mach  number,  typically  about  5  or  6,  the  entire 


isomagnetic  jump  disappears  and  the  shock  structure  is  taken  over  by  the 
magnetic  pedestal.  Hence,  we  feel  that  a  key  element  in  a  fluid 
description  of  shocks  in  this  intermediate  range  is  the  breakdown  of  quasi¬ 
neutrality.  Once  quasi -neutrality  is  violated,  an  additional  fluid 
dissipation  mechanism,  longitudinal  resistivity  becomes  possible.  In  the 
next  two  sections  we  examine  the  consequence  of  longitudinal  resistivity. 
Although  the  actual  dissipation  mechanism  is  considerably  more  complicated 
and  non  fluid  like  (see  Sec.  IV),  a  longitudinal  resistivity  model  can 
explain  many  features  of  the  experimental  results.  Specifically  it  hinges 
on  the  breakdown  of  quasi-neutrality;  it  provides  for  the  proper  electron 
and/or  ion  heating;  it  predicts  an  isomagnetic  potential  jump  on  the  crest 
of  the  shock  wave;  and  it  gives  the  shock  structure  only  in  an  intermediate 
range  of  Mach  numbers  between  the  critical  Mach  number  for  transverse 
resistivity  and  the  onset  of  strong  ion  reflection.  Thus,  longitudinal 
resistivity  is  probably  as  good  an  approximation  to  the  dissipation 
mechanism  that  one  can  obtain  within  the  limits  of  a  conventional  fluid 


description 


III.  Shocks  with  Transverse  and  Longitudinal  Resistivity 


In  this  Section,  we  discuss  shocks  formed  by  longitudinal  resistivity. 
Before  doing  so,  however,  in  Section  III. A  we  set  the  stage  by  reviewing 
the  structure  of  shocks  with  transverse  resistivity.  Section  III.B 
discusses  the  effect  of  longitudinal  resistivity,  and  Section  III.C  sets 
out  the  procedure  for  calculating  the  shock  structure  in  weakly  magnetized 
plasmas . 
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..  Review  of  Shocks  Formed  by  Xransver3e  Resistivity 


To  start  we  formulate  the  equation  for  quasi-neutral  transverse 
resistive  shocks.  We  neglect  both  transverse  and  longitudinal  electron 
inertia,  and  the  motion  of  the  ions  in  the  transverse  direction.  Assuming 
that  the  shock  propagates  in  the  x  direction,  the  magnetic  field  is  in  the 
z  direction,  and  upstream  parameters  have  a  zero  subscript,  then  the 
electron  and  ion  density  equation  gives  the  result 


nV  -  n  V  , 
o  o’ 


where  a  velocity  or  electric  field  without  a  direction  subscript  means  a 
velocity  or  electric  field  in  the  x  direction.  Then  the  electron  momentum 
equation  in  the  x  and  y  direction  give  respectively  equations  for  E  and 
Vey.  They  are 


--j-n  T  -  B 
ne  dx  e  c 


- 2.  (b - —  B] , 

mvc  '  o  n  1  ’ 


where  m  is  the  electron  mass,  v  the  effective  collision  frequency,  and  the 

resistive  force  density  on  the  electrons  is  -  nmvfv  -  V,  1  =  -  nmv  V 

v  ey  iy  •  ey 

since  we  assume  V  «  0  (unmagnetized  ions) .  Also  we  have  used  the  fact 
V  y 

that  Ey  -  V  “  “  V 

To  write  out  the  electron  temperature  equation  we  first  assume  that 
the  electrons  have  y  «  5/3.  Then  we  must  also  specify  what  fraction  of  the 


resistive  heating  goes  into  the  electrons  and  what  fraction  goes  into  the 
ions.  In  terms  of  the  ion  y,  the  total  energy  equation  for  electrons  plus 
ions  is 


-  neV  E  , 
ey  y» 


whereas  the  total  energy  equations  for  the  electrons  and  ions  separately 


j-  j  nVT 

dx  2  e 


-  neV  »E  +  A 


(1  nMV3  +  Jfj.  nVTj  -  neVE  -  A, 


where  M  is  the  ion  mass.  The  term  A  denotes  an  energy  exchange  between 
electrons  and  ions,  and  its  value  depends  on  what  fraction  of  the  resistive 
heating  goes  to  electrons  and  what  fraction  goes  to  ions.  To  determine 
this,  one  must  write  out  the  temperature  equations  for  electrons  and 
ions.  To  do  so,  take  the  dot  product  of  the  electron  momentum  equation 
with  V_  and  subtract  the  result  from  the  electron  energy  equation  to  get 


~  nVT  -  V  nT  -  A  +  nmv\ 
2  dx  e  dx  e 


Then  multiply  the  x  component  of  the  ion  momentum  equation 


by  V,  subtract  the  results  from  the  ion  energy  equation  to  obtain 


If  a 
then 


4  ~  nVT  -  V  ~  nT.  -  -  A. 
2  dx  i  dx  i 


(8) 


Is  the  ratio  of  ion  resistive  heating  to  electron  resistive  heating. 


or 


a 


A 


A  +  nmvV^ 
ey 


(9a) 


A 


a 

1  +  a 


nmv 


ey 


(9b) 


The  quantity  a  depends  on  the  physical  process  causing  the  resistive 
heating.  Classical  resistivity  gives  a  ~  m/M.  Anomalous  resistivity  gives 
an  a  which  depends  upon  the  particular  mechanism.  Generally  a  «  1;  for 
Instance  if  the  resistivity  is  caused  by  an  ion  acoustic  instability, 
a  -  0.1-0. 2.  Finally,  these  equations  are  closed  by  Ampere's  law 


dB  m  4irne  v 
dx  c  ey 


(10) 


In  order  to  see  where  the  critical  Mach  number  is,  a  standard 
procedure  is  to  linearize  Eqs.  (1-10)  about  the  upstream  or  downstream 
state.  Assuming  that  ■  k,  and  an  unsubscrlpted  variable  denotes  either 
the  upstream  or  downstream  state,  there  are  only  two  roots  for  k, 


k  -  0 


(11a) 


2  ,.2  2  yTi  5  Te 

k  -  “?el  f  ~  A  ~  M  "7M, 

-  v2  ill  ’ 

V  M  3  M 


(11b) 


where  V is  the  Alfven  speed.  In  our  convention,  V  >  0  so  that  a  state 


with  k  >  0  diverges  from  the  upstream  or  downstream  state.  Notice  that  k 


is  always  independent  of  a*  The  reason  is  that  all  a  terms  are  multiplied 


by  V  which  vanish  on  linearization.  Further,  notice  that  there  are  three 
ey 


ranges  of  velocity.  First,  if 


2  2  ^T.  5  ^e 

V  >  V  +  — +  —  — 
'  A  M  3  M 


k  >  0  so  the  linearized  solution  diverges  away  from  the  state.  Equation 


(12)  corresponds  to  the  upstream  state  of  the  shock  wave  for  which  the  flow 


speed  must  exceed  the  magnetosonlc  velocity.  The  second  regime  corresponds 


^i  5  ^e  2  2 

+  .d  <  v  <V  + 
1  3  M  s  A 


*t.  c  t 

1  1  ,  5  _e 

M  3  M  ’ 


where  k  <  0.  Since  the  flow  speed  is  less  than  the  magnetosonlc  speed,  it 


corresponds  to  the  downstream  state,  and  since  k  <  0,  perturbations 


converge  into  the  downstream  state.  This  then  is  a  standard  (transverse) 


resistive  shock  having  a  subcritical  Mach  number.  The  final  regime 


corresponds  to 


having  k  >  0.  Since  perturbations  now  diverge  away  from  the  downstream 
state,  there  can  be  no  shock  in  the  regime  defined  by  Eq.  (14)  if 
transverse  resistivity  is  the  only  dissipation  mechanism.  This  then 
defines  the  critical  point  above  which  no  shock  can  form. 

To  continue,  let  us  write  out  the  equations  for  a  quasi-neutral  shock 
as  a  system  of  nonlinear  first  order  differential  equations,  as  is 
appropriate  for  standard  numerical  integration. 

The  equations  for  B  and  Vgy  are  already  expressed  in  the  proper  format 
in  Eqs.  (10)  and  (3).  Equations  (6a)  and  (6b)  can  easily  be  reduced  to 


„  nmv 
2  r _ ey 

q  L  n  V 


\r  +  a  t 

ey _  _e  dn-, 

n  V  n  dx '  * 

o  o 


(16) 


and 


dT.  i  *  T.  . 

_ i  m  1  r  A  __i  dn-| 

dx  y-1  *-  n  V  n  dx ' 

o  o 


(15) 


but  now  dn/dx  appears  on  the  right  hand  side.  To  solve  for  dn/dx,  use  the 

ion  momentum  equation,  insert  the  electric  field  from  Eq.  (2),  use  Eqs. 

dTE  dTi  dn 

(15)  and  (16)  for  and  ,  and  solve  for  ^  . 


2  2 

V  ,  n  mvV 

ey  R  2  _ ey 

ne  B  >• 

c  3  n  V 


The  result  is 


This  can  be  used  on  the  right  hand  side  of  Eqs.  (15)  and  (Id).  Notice 
.  dT.  dT  , 

that  and  E  all  become  infinite  at  V  -  +  -j  Tg)/M.  Thus, 

2 

quasi-neutral  shocks  formed  by  transverse  resistivity  only  exist  if  V  < 

+  y  T^j/M  everywhere.  If  the  shock  is  monotonic  and  y  *  y  (or  more 
generally,  if  yg  *  Y^)»  this  critical  Mach  number  can  be  determined  by  the 
Rankine-Hugoniot  conditions  only.  However,  if  the  shock  is  not  monotonic, 
and/or  yg  *  Y^»  the  critical  Mach  number  depends  on  the  details  of  the 
oscillatory  structure  and  how  the  resistive  heating  is  partitioned  between 
electrons  and  ions. 

A  key  question  then,  for  the  description  of  shocks  in  the  intermediate 

regime,  is  the  physical  significance  of  the  singularity.  One  might  think 

that  the  singularity  indicates  the  onset  of  ion  reflection,  but  this  is  not 

the  case.  To  see  this,  imagine  that  T^  <<  Tg.  Then  ion  reflection  occurs 
1/2 

where  V  -  (T.  /M)  »  0.  However,  the  singularity  occurs  at  a  much  higher 

5Te  1/2 

velocity,  V  •  (y^— )  .  The  significance  of  the  singularity  lies  not  in 

ion  reflection,  but  rather  in  the  breakdown  of  the  quasi-neutral 

dE 

assumption.  This  is  clear  from  the  fact  that  both  E  and  yy  ,  as  expressed 
in  Eq.  (3),  become  infinite  at  the  singularity,  implying  a  breakdown  of 
quasi-neutrality.  Thus,  in  the  intermediate  regime,  charge  separation 


begins  to  play  an  Important  role  in  the  shock  structure. 


Shocks  with  Longitudinal  Resistivit 


In  this  subsection  and  the  next,  we  discuss  how  shocks  are  formed  by 
longitudinal  resistivity.  If  quasi-neutrality  is  no  longer  assumed,  the 
steady  state  ion  fluid  equations  are 


n.V.  -  n  V 
i  i  oo 


(18) 


a?  (V^i  +  niTi)  "  nieE  “  ni°Vvi 


V’ 


(19) 


and 


It  (t  "imi  +  A  "iviTi) '  “iVE  +  if;  "i"uv2 


ey 


+  TT^V“tfVVe>2- 

l 


(20) 


where  a#  and  v  are  defined  analogously  to  a  and  v.  Notice  that  we  allow 

X  X 


different  values  for  longitudinal  and  transverse  dissipation  and  energy 
partition.  This  is  reasonable  because,  as  we  will  see,  these  two  types  of 
dissipation  arise  from  completely  different  physical  processes. 

The  fluid  equations  for  the  electrons  are 


n  V  -  n  V  , 
e  e  oo 


(21) 


-r—  n  T  - 
dx  e  e 


V  B 
ey 


-  n  eE  -  n  e  — —  +  n.mv.fv.  -  V  ), 
e  e  c  ii'i  e' 


(22) 
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U  J 
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n  V  T 
e  e  e 


-  n  V  eE  -  n  V  eE  +  — 

e  e  e  ey  y  1  + 


—  n.mvV 
a  i  ey 


1  + 


ninvi(vi '  VJ: 


and 


V 

ey 


eV 

(B 

mvc  o 


These  are  supplemented  by  Maxwell's  equations 


dB 

dx 


4irn 

- ®  V 

c  ey 


and 


dE 

dx 


(23) 


(24) 


(25) 


(26) 


As  in  the  previous  section,  we  begin  by  linearizing  these  equations 

about  the  upstream  or  downstream  state.  As  before,  there ^is  no  resistive 

2  ^  n , 

heating  in  linear  theory,  so  f  ■  -j  T  —  and  *  (y-l)^  —  .  Of  course, 

o 

»  -  V  n  /n  and  m  -  V  n./n.  Then  linearizing,  a  straightforward 
e  e  1  1 


calculation  gives  the  result 


separates.  The  low  k  root  is  still  given  by  Eq.  (lib).  We  denote  this  the 


Magnetic  Mode.  To  get  a  simple  form  for  the  high  k  roots,  we  further 
2  2 

assume  that  moj  /  Mlc_v  V  «  1,  and  u>  /v  <<  1.  Since  the  length  of  a 

c  u  o  c  i 

resistive  shock*-^  is  typically  10c/u>pe,  an  approximate  order  of  magnitude 

£  _  w  *.k«  AAlT4A4An  £  4  r.  ..  ..  in..  17  y  ^  TU«  C4 


v  V  _  » 

-T-  (y  Te  -  -  MV 

u 

T-  (^1  -  mo)  1 


TT  (''I1-flv4)  "u; 


1  Trt(^[4Ie+’'Ii-MVol)2  + 


+  - 11,2 )  WTi  -  MVo)'1/2 


(28) 


where  ug  is  the  electron  thermal  velocity  ug  -  (Te/m)^^. 

Notice  that  the  wave  number  given  by  Eq.  (29)  is  of  order  kQ.  Also, 
the  magnetic  field  completely  decouples  from  this  large  k  mode,  so  it  is 
entirely  electrostatic  in  nature.  Depending  on  the  values  of  velocity  and 
temperature,  the  linearized  modes  structure  themselves  as  shown  in  Table  I. 

Now  consider  how  to  integrate  Eqs.  (18)— (26)  to  solve  for  the  shock 
structures.  Let  us  say  we  start  upstream  and  integrate  toward  the 
downstream  state.  Referring  to  Table  I  the  upstream  plasma  is  in  Class  A, 
since  the  flow  speed  must  exceed  the  magnetosonic  speed.  Hence  there  are 
two  linearized  perturbations  which  increase  away  from  the  upstream  state. 
The  relative  amplitude  between  the  two  is  arbitrary  at  this  point.  Now  let 
us  assume  the  downstream  plasma  is  in  Class  B  or  C.  In  each  case  there  are 
two  linearized  eigenfunctions  which  converge  into  the  downstream  state,  and 
one  which  diverges.  Thus  the  relative  amplitude  of  the  two  linearized 
perturbations  upstream  must  be  chosen  so  as  to  eliminate  the  divergent 
perturbation  downstream.  This  is  then  a  nonlinear  eigenvalue  problem,  and 
this  concept  has  been  applied  to  both  conventional  shocks^  and  reflected 
ion  shocks. ^ 
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Table  I 


Class 


Velocity  Regimes  and  linearized  modes  in  each  regime. 

Velocity  Range _ Modes _ 


MV?  +  t  T  +  vT,  <  MV2 
A  3  e  i 


One  electrostatic  mode  and  the 
magnetic  mode  diverge  from  the 
state  as  x  increases  in  the 
direction  of  flow.  The  other 
electrostatic  mode  converges. 


y  Te  +  yTi  <  MV2  < 


ml  +  T  Te  +  *Ti 


Magnetic  mode  converges,  one 
electrostatic  mode  converges,  one 
diverges . 


YTt  <  MV2  <  y  Te  +  3Ti 


Magnetic  mode  diverges,  both 
electrostatic  modes  converge. 


MV  <  yT. 


Magnetic  mode  and  one  electric  mode 
diverge,  the  other  electrostatic 
mode  converges. 


k*  .v. ^ ;.vy- 


If  the  downstream  stats  is  Class  B,  the  scheme  is  trivial.  The 
divergent  downstream  perturbation  is  electrostatic  so  the  solution  is  just 
to  eliminate  electrostatics  and  calculate  shock  structure  based  on  quasi¬ 
neutrality.  If  the  downstream  state  is  Class  C,  the  problem  is  more 
complex.  Here  the  divergent  state  is  magnetic  and  both  convergent  states 
are  electrostatic,  so  there  is  necessarily  a  coupling  between  the  magnetic 
and  electric  structure  somewhere  in  the  shock.  In  the  next  subsection  we 
will  show  how  to  calculate  this  coupled  structure  in  a  simple  way  by 
exploiting  the  fact  that  the  length  scale  for  the  electric  structure  is 
very  small  compared  to  the  length  scale  for  magnetic  structure. 

Finally,  if  the  downstream  state  is  Class  D,  there  are  two  divergent 
states  downstream,  but  only  one  free  parameter  upstream,  so  that 
longitudinal  resistivity  does  not  allow  a  shock  to  form  in  this  case.  In 
the  next  subsection  we  will  discuss  the  meaning  of  the  singular  behavior  in 
this  regime. 
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Calculating  the  Doable  Structured  Shock 


As  we  have  seen  in  the  previous  two  subsections,  when  the  plasma  is 

making  the  transition  from  Class  B  to  Class  C,  the  electrostatic 

approximation  breaks  down  and  large  electric  fields  and  charge  separations 
are  set  up  by  the  flow.  If  the  plasma  were  in  Class  C,  and  there  were  no 
dissipation,  this  charge  separation  results  in  an  undamped  ion  acoustic 
wave.  If  it  were  in  Class  B  with  no  dissipation  the  charge  separation 

manifests  itself  as  an  ion  acoustic  soliton.  With  dissipation  in  the  form 
of  longitudinal  resistivity,  an  ion  acoustic  shock  can  form  so  that  the 

plasma  can  make  a  transition  from  Class  B  to  Class  C.  Thus  near  this 
point,  a  narrow  ion  acoustic  shock  is  embedded  in  the  broader  magnetic 
structure.  This  then  corresponds  to  the  isomagnetic  potential  jump 
discussed  in  Ref.  23. 

The  choice  of  the  electrostatic  wave  amplitude  upstream  is  then 

equivalent  to  a  choice  of  the  position  of  the  ion  acoustic  subshock  within 

the  broader  magnetic  structure.  As  we  have  seen,  for  a  Class  C  plasma 

there  is  no  magnetic  perturbation  which  converges  into  the  downstream  side 

2  2 

of  the  ion  acoustic  subshock.  This  means  for  the  case  that  kn  >>  oj  V/vc  , 

u  pe 

where  the  magnetic  structure  and  ion  acoustic  wave  decouple,  the  ion 
acoustic  wave  must  be  at  the  position  where  the  magnetic  field  has  its 
downstream  value. 

Hence  the  following  scheme  for  computing  the  shock  structure  suggests 
itself.  First  integrate  the  equations  for  a  quasi-neutral  resistive  shock, 
as  written  in  Section  III. A  until  one  reaches  the  position  where  the 
magnetic  field  has  the  downstream  value  predicted  by  the  Rankine-Hugoniot 
conditions.  We  find  that  at  this  point,  the  densities  and  temperatures  are 
always  below  the  downstream  values,  and  also  MV  >  5/3  Tg  +  yT^.  At  this 
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point  simply  set  the  magnetic  field  equal  to  its  constant  value,  set  Vgy  • 
0  and  solve  for  the  ion  acoustic  shock  in  the  unmagnetized  plasma.  The 
equations  of  the  ion-acoustic  shock  are 


s  (r  “X 


n.V  eE  +  -r-r- —  n.mv.fV, 
i  i  1+a^  i  JA  i 


(29) 


The  ion  density  and  momentum  equation  are  as  in  Eqs.  (18)  and  (19).  For 
the  electrons 


n  V 
e  e 


n  V  , 
o  o 


(30) 


■rr—  n  T  = 
dx  e  e 


-  n  eE  +  n.mv  fV.  -  V  V 
e  i  a  ^  i  e 1 


(31) 


4 —  4-  n  V  T  - 
dx  2  e  e  e 


-  n  V  E  +  rr^-r -  n.mv.fV.  -  V2  ] , 

ee  1+a  l  i  eJ 


(32) 


and  Poisson's  Equation  is  still  as  given  in  Eq.  (26).  At  the  onset  of  the 
ion  acoustic  shock,  the  electric  field  can  be  initialized  by  taking  its 
value  from  the  quasi-neutral  solution  at  the  switch  point.  Of  course,  the 
exact  value  that  E  is  initialized  at  is  not  very  important  because  E 
reaches  much  higher  values  in  the  center  of  the  transition  region. 

Then  Eqs.  (18,  19,  26,  29-32)  are  integrated  from  the  start  of  the  ion 
acoustic  shock  on  through  to  the  downstream  state.  It  is  easy  to  show  that 
this  procedure  preserves  the  overall  Rankine-Hugoniot  conditions;  this  is 
shown  in  Appendix  A. 

To  proceed,  let  us  examine  the  singular  nature  of  these  equations  by 

dTi 

writing  out  equations  for  V^.  Eliminating  from  the  energy  equation 


(Eq.  (29))  and  inserting  it  in  the  momentum  equation,  Eq.  (19),  we  find 


(33) 


dVt  eViE  -  “\(V1  -  VeH^t  -  ^  TT^  (vi  -  ven 

d*  ;MVE  -  vTj 

Now  the  only  singularity  is  at  *  yT^/M.  It  is  especially  easy  to  draw 

insight  from  this  for  the  case  y  ■  3,  one  dimensional  ions.  As  shown  in 

Ref.  10,  the  fluid  equations  are  those  derived  from  a  waterbag  ion 

1/2 

distribution  for  which  the  thermal  width  is  (31^/M)  .  Then  the 

2 

singularity  at  V  *  31^/M  simply  means  that  an  ion  at  the  outer  edge  of  the 

waterbag  distribution  function  has  been  reflected  (its  velocity  has  gone  to 

zero  in  the  shock  frame).  The  obvious  generalization  then  is  to  regard  the 
2 

singularity  at  V  *  yT^/M  as  the  generalization  of  the  condition  for  ion 
reflection  to  two  or  three  dimensional  ion  distributions.  Hence  this 
singularity  does  have  a  simple  interpretation  in  terms  of  reflected  ions. 
In  Ref.  10,  it  was  shown  that  for  ion  acoustic  shocks,  one  could  add  a 
reflected  fluid  upstream.  The  density  of  this  reflected  fluid  then  becomes 
an  additional  upstream  parameter,  which  can  be  varied  so  as  to  eliminate 
another  divergent  downstream  solution.  This  might  also  be  possible  to  do 
for  the  magnetosonic  shocks  which  we  consider,  but  it  will  not  be  examined 
further  here. 

To  summarize,  we  have  shown  that  longitudinal  resistivity  can  give 
rise  to  double  structured  shocks  in  an  intermediate  range  of  Mach  number. 
These  are  where  the  downstream  state  is  Class  C,  as  defined  by  Table  1.  At 
the  upper  limit  of  this  range,  the  shock  structure  once  again  becomes 
singular,  but  this  singularity  can  now  be  interpreted  in  terms  of  reflected 
ions. 
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In  the  previous  section,  we  discussed  how  longitudinal  resistivity 


forms  the  transverse  shock  in  the  intermediate  Mach  number  regime  (Class 


C).  The  idea  is  that  as  the  Mach  number  increases,  large  charge  separation 

is  generated,  so  an  ion  acoustic  wave  is  set  up.  As  this  wave  damps  by 

longitudinal  resistivity,  the  shock  forms,  as  discussed  in  Section  II. 

However,  longitudinal  resistivity  is  only  a  fluid  model  for  the 

dissipation,  the  actual  dissipation  mechanism  is  electron  or  ion  Landau 

damping,  which  is  inherently  non-fluid  like.  Nevertheless,  it  is  an 

extremely  powerful  dissipative  mechanism.  If  the  ions  are  cold  the  damping 

distance  due  to  electron  Landau  damping  is  extremely  small,  of  order 

(M/m)*^  k  This  implies  a  v  of  order  oj  .  If  the  ions  are  warm,  the 
D  i  pe 

Landau  damping  rate  is  larger,  implying  still  larger  values  of 

As  the  longitudinal  resistivity  appears  first  in  the  momentum 
equation,  we  proceed  by  calculating  the  term  in  the  electron  momentum 
equation  which  arises  from  electron  Landau  damping.  This  means  getting  the 
correction  to  the  pressure  tensor.  To  do  so  we  assume  the  plasma  is 
infinite  and  homogeneous  work  in  the  (u>,  k)  domain,  and  then  write  a 
linearized  fluid  equation.  The  only  term  which  is  not  completely 

straightforward  is  the  pressure  tensor.  As  usual,  we  denote  perturbed 
quantities  with  a  superscript  tilda.  Then 

!ie 

g  eE  3v 

e  m  -iui  +  ikv 
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(34) 


and 


eE  3v 
M  -ioj  +  ikv  ’ 


where  now  a  v  denotes  a  particle  velocity  in  a  Vlasov  model  and,  as  usual, 
a  capital  V  is  the  fluid  velocity.  The  total  perturbed  electron  density  is 


-  eE  ,  3v 

ne  =  m  J  -ico  +  kv 


ikm  u 


[1  +  f  -  f 

L  1  ~u)  +  kv  e 


_  2  2 

where  we  have  assumed  Maxwellian  electrons,  fg  ■  (/2it  u^)  exp-(v  /2ue)  and 
we  work  now  in  the  plasma  rest  frame.  Taking  the  resonant  part  of  the 
second  term  in  the  brackets  of  Eq.  (36),  that  is 


Im  -u  VkV  ' 


we  find 


~  eE  r-  /ir-vl/2  u>  ■) 

e  77 — ?  w  TCTTr' 

ikm  u  e 
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To  continue,  we  would  like  to  find  the  perturbed  part  of  the  total 
electron  momentum  flux 


Pe  -  m  /  dv  v2  Fe  -  -  *  /  dv  v3  _r-  k~  •  (39) 

u 

e 

The  resonant  part  Is  small  by  a  factor  of  (j~~ )  and  we  neglect  it,  so  that 


p  .  _  eE  n 

e  “IF"  * 
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To  obtain  an  equation  of  state,  we  need  the  perturbed  pressure  in  terms  of 
the  perturbed  density.  From  Eqs.  (38)  and  (40)  we  find 


(41) 


The  force  on  the  electrons  is  -ikt^.  The  first  term  in  the  square  brackets 
simply  gives  -  T  —  n  ,  the  pressure  gradient  force  on  an  isothermal 
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plasma.  The  second  term,  which  is  out  of  phase  by  ir/2  is  the  additional 
dissipative  contribution  to  the  electron  force.  Using  the  fact  that  the 
linearized  density  equation  gives  the  result  n if  -  m/k  ne,  we  find  the 
additional  force  density  on  the  electrons  is 
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Modeling  this  additional  retarding  force  as  a  longitudinal  resistivity  in 
the  subshock,  we  find 


v 


l 


u> 


pe 


(43) 


where  we  have  assumed  that  within  the  subshock  k  ~  k^. 

Although  the  dissipation  from  Landau  damping  is  quite  simple  in  the  k 
domain,  it  is  very  complicated  and  nonlocal  in  the  x  domain.  In  the  k 
domain,  the  dissipative  part  of  the  perturbed  pressure  is  the  product  of 
-i  nm  ue  and  In  the  x  domain,  it  is  the  convolution  of  the 

Fourier  transforms  of  the  two  factors.  Thus  the  dissipative  part  of  the 
perturbed  pressure  in  the  x  domain  is 


28 


",  Mr*  -Tm  Mm  ^m  ”  •  w  •  “  *  .  '  • 


V  (x') 

n  m  u  t/^r  P  /  dx'  — — 
e  2  J  x-x 
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where  P  denotes  the  principal  value.  Thus  the  actual  dissipation  is 
nonlocal  and  quite  complicated. 

Now  we  consider  the  dissipative  force  on  the  ions.  An  exact  analog  of 
the  calculation  leading  up  to  Eq.  (42)  gives  the  result  that  the 
dissipative  force  density  on  the  ion  fluid  is  given  by 


*My  TkT^r'ui/kV 


(45) 
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For  the  case  of  positive  phase  velocity,  |  <  0  for  a  thermal  ion 

distribution,  so  that  the  force  is  still  a  retarding  force.  Assuming 
id/ k  is  roughly  constant  as  a  function  of  k,  the  perturbed  ion  pressure  in 
the  x  domain  has  the  same  form  as  given  in  Eq.  (45).  Depending  on  the  ion 
temperature,  the  corresponding  to  the  ion  Landau  damping  can  be  large  or 
small.  If  the  electron  temperature  is  comparable  to  the  ion  temperature, 
the  damping  length  is  comparable  to  k^,  implying  that  <  (— )  “pe’ 

Hence  for  either  electron  or  ion  dissipation,  values  for  are  quite 
large,  easily  large  enough  that  the  electrostatic  potential  is  highly 
localized  in  the  overall  magnetic  structure.  Furthermore,  the  dissipation 
is  not  anomalous,  in  the  sense  that  it  does  not  arise  from  any  instability, 
as  does  the  dissipation  which  generates  the  transverse  v. 

Since  ion  Landau  damping  heats  the  ions,  and  electron  Landau  damping 


heats  the  electrons,  a  reasonable  guess  for  a.  is 


/2rr  MV 
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where  we  have  taken  aj/k  «  V,  the  flow  speed. 

We  now  continue  this  section  by  discussing  other  aspects  of  the  Landau 
damping  model.  As  complicated  as  say  Eq.  (44)  for  Pgd  is,  it  is  still  a 
fairly  coarse  approximation  to  the  actual  nonlocal  dissipation  for  several 
reasons.  First  of  all,  Eq.  (44)  is  derived  from  linear  theory,  whereas  the  l 

shock  wave  is  inherently  nonlinear.  Secondly,  since  it  is  nonlocal, 

perturbed  velocities  everywhere  contribute  to  the  perturbed  dissipative 
pressure  at  each  point.  However,  only  a  portion  of  the  plasma  is 
electrostatic;  a  large  part  of  it  is  dominated  by  the  magnetic  structure, 
and  this  was  not  treated  correctly  at  all  in  the  calculation  of  Pgd. 

Third,  the  Fourier  transform  technique  used  to  calculate  Pg  works  only 
if  Pg  is  either  periodic  or  square  integrable.  However,  a  shock  profile  is 
neither  so  it  is  not  clear  just  how  valid  the  Fourier  representation  is. 

For  instance,  if  Vg(x)  in  Eq.  (44)  has  a  typical  shock  structure,  then  the 
integral  in  Eq.  (44)  diverges  for  large  x.  Thus  the  linear  theory  in  the 
Fourier  domain  gives  an  infinite  perturbed  pressure  at  large  x,  that  is, 
the  region  of  x  furthest  from  the  shock.  It  is  the  divergence  which  is  at 
the  root  of  the  claim  by  Pfirsch  and  Sudan  that  Landau  damping  cannot 
form  a  steady  state  shock.  We  feel  that  this  divergence  is  clearly 
unphysical,  so  either  nonlinear  contributions  play  an  important  role  in 
determining  Pg  or  else  the  Fourier  transform  scheme  must  be  modified. 

Hence,  since  the  actual  dissipation  mechanism  is  so  complicated,  a 
simple  fluid  approximation  to  it  could  be  extremely  useful.  The  only  four 
possible  conventional  dissipation  mechanisms  are  thermal  conductivity, 
viscosity,  transverse  resistivity,  and  longitudinal  resistivity.  Thermal 
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conduction  can  almost  certainly  be  eliminated  as  a  mechanism  to  form  the 
shock  because  there  does  not  seem  to  be  any  way  that  diffusing  temperature 
can  slow  the  forward  motion  and  turn  the  upstream  kinetic  energy  into 
downstream  temperature.  Viscosity  can  form  the  shock.  However,  there  is 
no  maximum  Mach  number  for  a  viscous  shock,  no  way  that  It  Is  related  to 
the  breakdown  of  quas  1-neutrality ,  and  no  way  that  it  is  related  to  ion 
reflection.  Since  all  of  these  are  known  to  be  important  in  experiments, 
viscosity  does  not  appear  to  be  the  best  fluid  description  of  the 
dissipation  either.  Transverse  resistivity  is  not  only  a  good 
approximation,  in  all  likelihood  it  is  the  dissipation  mechanism,  since 
transverse  electron-ion  streaming  can  give  rise  to  instabilities  which 
exchange  momentum  between  electrons  and  ions.^  However,  as  we  have  seen, 
transverse  resistivity  can  only  form  the  shock  for  Class  B  plasmas. 

This  leaves  longitudinal  resistivity.  While  the  actual  dissipation 
mechanism  is  considerably  more  complex  than  this,  it  is  a  very  good 
approximation  to  the  physics.  It  is  directly  related  to  the  breakdown  of 
quasi-neutrality,  it  has  an  upper  Mach  number  limit  corresponding  to  ion 
reflection,  it  predicts  a  very  narrow  isomagnetic  potential  jump  at  the 
shock  crest,  and  it  has  a  lower  Mach  number  limit  corresponding  to  the 
maximum  Mach  number  for  transverse  resistivity.  All  of  these  features  are 
confirmed  by  the  laboratory  experiments  we  have  cited.  Furthermore,  while 
longitudinal  resistivity  is  simple,  the  two  parameters  which  characterize 
it,  v.  and  a.  can  be  estimated  from  a  coarse  approximation  to  the  actual 
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dissipation  mechanism. 

Let  us  conclude  this  section  by  discussing  an  important  additional 
constraint  on  the  physics.  As  we  have  seen,  our  basic  model  is  that  the 
flow  steepens  up  until  an  ion  acoustic  wave  is  generated.  The  Landau 


damping  of  this  ion  acoustic  wave  then  provides  the  dissipation  which  forms 
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the  shock.  However,  for  the  theory  of  Landau  damping  to  be  valid,  the 
trapping  time  (for  say  the  ions) 

tt  ~  k  <47> 


where  $  is  the  electrostatic  potential,  must  be  longer  than  the  correlation 
time  for  a  particle  with  velocity  v. 


Z  |E(k)|25(kv-u>) 
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Z  |E(k) | 2 
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(48) 


However,  if  the  shock  is  simply  a  steady  state  potential  profile,  t  =■  00 

c 

for  particles  at  the  shock  velocity.  Hence  the  entire  concept  of  Landau 
damping  breaks  down  and  the  ion  dynamics  is  dominated  by  trapping  in  a 
stationary  potential.  Thus  in  order  for  the  Landau  damping  model  to  be 
valid,  either  the  shock  must  have  a  temporally  fluctuating  structure,  or 
else  the  total  time  of  the  experiment  must  be  less  than  a  trapping  time. 
It  is  very  likely  that  the  shocks  do  have  a  fluctuating  structure  if  only 
because  in  the  broader  part,  dominated  by  transverse  resistivity,  the 
resistivity  is  probably  caused  by  Debye  length  scale  streaming 
instabilities.  The  spectrum  of  unstable  waves  will  almost  certainly 
overlap  the  electrostatic  subshock  and  give  rise  to  overall  fluctuations  in 


V.  Results 


In  this  section  we  present  results  of  the  calculations  outlined  in 
Section  III.  There  are  two  classes  of  calculations.  First,  we  consider 
parameters  like  those  examined  experimentally  in  Ref.  22.  Here  we  model 
the  longitudinal  resistive  heating  as  going  entirely  into  the  ions,  as 
would  be  the  case  for  strong  ion  Landau  damping  (that  is  not  much  less 
than  T0).  Secondly,  we  consider  a  nitrogen  plasma  with  parameters  like 
that  given  in  Ref.  21.  More  important,  however,  we  consider  no  resistive 
ion  heating. 

A.  The  Case  of  Strong  Ion  Heating 

In  this  case,  we  consider  a  shock  in  a  hydrogen  plasma  with  parameters 

as  given  in  Ref.  22.  That  is,  the  upstream  density  is  4.6x10^  cm-^  and 

the  upstream  field  is  730G.  The  ion  y  is  taken  as  equal  to  5/3,  and  Tg  » 

T^  -  1  ev.  As  is  conventional  in  resistive  shocks  we  assume  a  length  of 

roughly  10  c/co  so  that  cv/ui  V.  *  10  where  we  use  upstream  values  for  the 
pe  pe  A 

plasma  frequency  and  Alfven  speed.  Although  it  is  often  assumed  that  for 
transverse  resistive  shocks  there  is  no  resistive  ion  heating,  it  is  not 
easy  to  envision  a  transverse  streaming  instability  which  gives  none. 
Also,  a  small  amount  of  resistive  ion  heating  would  be  extremely  difficult 
to  measure,  and  it  might  in  fact  be  very  important  in  determining  the 
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resistive  shock  width.  In  addition,  the  ion  temperature  greatly  affects 
the  ion  Landau  damping  rate  in  the  electrostatic  portion  of  the  shock. 
Accordingly  we  take  a  *  0.2,  so  there  is  5  times  as  much  resistive  electron 
heating  as  resistive  ion  heating  in  the  transverse  part  of  the  shock. 


Since  the  transverse  resistive  shock  always  has  a  monatonic  structure, 
the  critical  Mach  number  for  transition  to  the  Class  C  downstream  state 
can  be  determined  from  the  Rankine-Hugoniot  conditions.  For  the 
parameters  cited,  this  critical  Mach  number  turns  out  to  be  2.64.  (Here 
the  Mach  number  is  defined  as  the  shock  velocity  divided  by  the  upstream 
magnetosonic  speed).  For  Mach  numbers  below  this,  the  shock  is  quasi¬ 
neutral  and  its  structure  is  determined  as  specified  in  Section  III. A. 

Above  this  Mach  number,  the  quasi-neutral  assumption  breaks  down.  For 
larger  Mach  numbers,  the  shock  structure  is  calculated  as  described  in 
Secs.  III.B  and  C.  We  assume,  as  seems  consistent  with  experimental  data, 
that  in  the  electrostatic  portion,  ion  Landau  damping  dominates.  Thus  we 
take  o.  *  <*»  and  v  (— .  The  upper  critical  Mach  number,  the 
transition  from  Class  C  to  Class  D  cannot  be  determined  from  Rankine- 
Hugoniot  conditions  alone  for  two  reasons.  First  of  all  the  electrostatic 
part  of  the  shock  generally  is  not  monatonic,  but  overshoots  the  final 
density,  and  then  decays  down  to  the  downstream  value.  Secondly,  even  if 
the  shock  were  monatonic,  the  Rankine-Hugoniot  condition  gives  only  Te  + 
T^,  whereas  the  transition  from  Class  C  to  Class  D  depends  only  on  T^. 
Therefore  this  upper  critical  Mach  number  depends  on  how  the  resistive 
heating  is  partitioned  between  electrons  and  ions.  For  the  parameters  we 
have  chosen,  we  find  that  this  upper  critical  Mach  number  turns  out  to  be 
about  4.46.  This  agrees  reasonably  well  with  the  upper  limit  Mach  number 
of  4.25  +  0.2  in  Ref.  22.  In  Figs,  la  and  lb  are  shown  magnetic  field  and 
density  profiles  for  a  Class  B  shock  with  m  ■  2.5.  Distances  are  given  in 
centimeters.  In  Fig.  2a  and  2b  are  shown  magnetic  field  and  density 
profiles  in  the  transverse  part  of  a  supercritical  shock  having  m  *  4.4. 


The  transverse  part  ends  when  the  magnetic  field  achieves  its  downstream 


value.  At  this  point  the  electrostatic  subshock  is  excited.  The  electron 
and  ion  density  as  a  function  of  x  in  this  region  is  shown  in  Fig.  2c. 
Again,  distances  are  in  centimeters.  In  all  calculations,  downstream 
values  are  carefully  compared  with  Rankine-Hugoniot  conditions  as  an 
independent  check  on  the  calculations. 

The  actual  ion  heating  downstream  was  measured  in  Ref.  22  by  measuring 
ne(Tg  +  T^)  magnetically,  and  then  measuring  ng  and  T0  by  laser 

scattering.  The  results,  along  with  our  calculated  results,  are  shown  in 
Fig.  3.  The  two  basically  agree,  although  the  experiment  shows  somewhat 
more  ion  heating  than  what  is  calculated. 

Another  important  experimental  result  from  Ref.  23  is  that  as  the  Mach 
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number  increases,  the  value  of  2e<})/M(Vo  -  V^)  decreases.  Here  <j>  is  the 
total  potential  jump  and  VQ  and  are  upstream  and  downstream  flow 

speeds.  This  goes  hand  in  hand  with  the  idea  that  as  the  Mach  number 
increases,  ions  are  slowed  down  not  only  by  the  electric  field,  but  also  by 
some  additional  dissipative  effect.  Our  dissipative  model,  longitudinal 
resistivity  also  provides  additional  ion  slowing  down.  In  Fig.  4  are  shown 
variations  of  this  parameter  with  Mach  number  as  measured  by  Eselevich  et 
al. ,  and  as  we  calculated  assuming  upstream  parameters  in  Ref.  22.  Again 
the  agreement  is  quite  good. 

B.  Shocks  with  No  Ion  Heating 

In  this  section  we  consider  shocks  with  no  ion  heating,  a  *  =  0. 

This  corresponds  to  zero  temperature  ions,  or  else  a  waterbag  ion 
distribution  in  one  dimension.  We  consider  the  latter,  so  y^  *  3.  Since 
the  electrons  and  ions  have  different  values  of  y,  the  Rankine-Hugoniot 
conditions  cannot  be  calculated  in  the  conventional  way.  However,  by 
making  use  of  the  fact  that  a  *  0,  and  the  ions  only  heat  adiabatically , 


the  conservation  relations  can  be  calculated.  This  is  done  in  Appendix 


C.  For  upstream  fluid  parameters,  we  use  those  in  Ref.  21,  a  singly 

ionized  nitrogen  plasma,  n  *  lO^"5,  B  =  2KG,  Tg  =  2  ev,  =  0.5  ev, 
cv/uipeVA  *  10,  and  =*  u>  ,  corresponding  to  electron-Landau  damping. 
Then,  Rankine-Hugoniot  conditions  show  that  the  Mach  number  for  transition 
from  Class  B  to  Class  C  is  M  =  2.68.  There  is  no  second  critical  Mach 
number  for  transition  from  Class  C  to  Class  D  predited  by  Rankine-Hugoniot 


simply  because  there  is  not  enough  ion  heating.  However,  the  shock  is  not 


monotonic,  but  oscillatory,  so  that  in  an  overshoot  it  is  possible  to  reach 
the  second  critical  point  V  ■  3T^/M.  We  find  this  singular  point  is 
reached  at  m  *  3.73.  In  Fig.  5a  and  5b  are  shown  magnetic  field  and 


density  profiles  versus  position  in  the  quasi  -meutral  part  for  a  shock 
with  m  m  3.5.  In  Fig.  5c  is  the  ion  density  versus  x  in  the  electrostatic 


VI. 


Conclusions 


Our  first  conclusion  is  that  the  experimental  evidence  strongly 
indicates  the  existence  of  a  significant  intermediate  range  of  Mach  numbers 
between  shocks  formed  by  transverse  resistivity  and  shocks  formed  by 
reflected  ions.  In  this  range,  according  to  Ref.  23,  the  breakdown  of 
quasi-neutrality  plays  an  important  role  and  an  electrostatic  subshock 
forms  on  the  crest  of  the  magnetic  structure.  Second,  we  have  shown  that 
the  singular  behavior  of  the  quasi-neutral  fluid  equations  is  not  related 
to  ion  reflection  but  rather  to  the  breakdown  of  quasi-neutrality.  In  the 
intermediate  regime  then,  an  ion  acoustic  wave  train  is  set  up.  We 
speculate  that  the  Landau  damping  of  this  acoustic  wave  is  the  dissipation 
mechanism  which  forms  the  shock.  Third,  we  have  investigated  longitudinal 
resistivity  as  a  simple  fluid  model  for  the  dissipation  mechanism  for  shock 
formation.  We  conclude  that  it  is  a  reasonable  model.  It  is  important 
only  between  the  upper  Mach  number  for  transverse  resistivity,  and  the 
lower  Mach  number  for  a  fluid  picture  of  ion  reflection.  It  also  is 
directly  related  to  the  breakdown  of  the  quasi-neutrality  assumption,  and 
it  predicts  an  electrostatic  subshock  at  the  crest  of  the  magnetic 
structure. 
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Appendix  A 


In  this  appendix,  we  show  that  the  scheme  for  determining  the  shock 
structure,  as  specified  in  Section  III.C  satisfy  the  overall  conservation 
relations.  Three  of  the  conservation  relations,  that  for  mass,  momentum 
and  energy  are  satisfied  at  each  point  in  the  shock.  The  fourth  relation 
uB  »  UQB0  is  only  satisfied  downstream. 

Consider  now  the  conservation  equations  for  momentum  and  energy 
flux.  Each  flux  is  the  sum  of  a  fluid  and  electromagnetic  component.  In 
the  quasi-neutral  region,  the  electrostatic  part  of  the  momentum  flux  is 
neglected,  but  the  magnetic  part  of  the  momentum  flux,  B^/8ir  and  the 
electromagnetic  part  of  the  energy  flux  BEq  are  each  included  in  the 
formulation.  In  the  quasi-neutral  region  then,  n,  V,  and  T  are  determined, 
in  part,  so  that  particle,  momentum  and  total  energy  flux  are  conserved. 
At  the  switch  point,  B  takes  on  its  downstream  value  and  remains 
constant.  Therefore  in  the  electrostatic  region,  the  magnetic  part  of  the 
momentum  flux  and  electromagnetic  part  of  the  energy  flux  are  constant  and 
have  the  proper  downstream  value. 

In  the  electrostatic  region,  the  formulation  leads  to  a  conservation 

law  for  total  energy  flux  and  a  conservation  law  for  total  fluid  momentum 

E2 

flux  plus  electrostatic  momentum  flux  -  .  Downstream,  of  course,  E  =*  0 
and  the  fluid  momentum  and  energy  fluxes  are  as  given  by  the  Rankine- 
Hugoniot  condition. 

Now  consider  the  fluid  energy  flux  as  one  moves  back  upstream  toward 
the  switch  point  in  the  electrostatic  region.  Since  it  is  constant,  it 
matches  smoothly  onto  that  in  the  quasi-neutral  regime.  Next  consider  the 
momentum  flux.  Since  the  formulation  for  the  electrostatic  region  contains 
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a  momentum  flux  -  -g^-  ,  and  the  quasi-neutral  formulation  does  not,  the 
momentum  fluxes  do  not  quite  natch  at  the  switch  point  if  they  have  the 

E2 

same  values  downstream.  However,  if  -5—  is  a  very  small  part  of  the 

ott 

momentum  flux  at  the  switch  point,  the  error  in  momentum  flux  is  very 

small.  Using  the  expression  for  E  at  the  switch  point,  Eq.  (2),  and 

considering  say  the  ratio  of  electrostatic  to  thermal  momentum  flux,  we 
E^  2  2  2 

find  -s—  / nT  ~  feu  v/vc  k_l  <<  1.  In  all  of  our  calculations  of  shock 
8  it  e  1  pe  D' 

structure,  we  have  checked  that  the  electrostatic  component  of  the  momentum 
flux  is  negligible  at  the  switch  point,  so  that  total  momentum  is  conserved 
throughout  the  shock.  Of  course,  while  the  electrostatic  momentum  flux  is 
very  small  at  the  beginning  and  end  of  the  electrostatic  shock,  it  can  be 
very  large  at  points  in  between. 
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Appendix  B 


In  this  appendix,  we  calculate  Rankine-Hugoniot  conditions  for  the 

case  of  Ye  m  y  >  *  3,  but  where  the  ions  are  only  heated  adiabatically 

so  that  Ti  -  Tio  (IN2-  Then  denoting  upsteam  quantities  with  a  subscript 

1  o 

zero  and  downstream  quantities  with  no  subscript,  the  conservation 
equations  for  momentum  and  energy  flux  reduce  to  two  equations  for  n  and 
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In  writing  Eqs.  (Bl)  and  (B2)  we  have  used  the  fact  that  nQV0  =  nV  and  VQB0 
«  VB.  Now  each  equation  is  linear  in  Te,  so  solving  for  Te  from  one  and 
inserting  in  the  other,  we  find 


MV 

MVo(1  +  W1  -  “  )  +  Ti2-  (l  -  “2)  +  -  “)Te0  - 

U 


(B3) 


n  2  2 

where  u  *  —  and  V.  ■  B  /4irn  M.  Factoring  out  the  obvious  root  u 
n  Ao  o  o 

o 

arrive  at  a  cubic  equation  for  u 


1,  we 


2  2 
MV  ,  MV* 

- -  +  T.  fl  +  u  +  u2)+  (1  +  u)  +  T  - 

u  io''  '  2  v  1  eo 

MV2 

!„[-»  (i  +  u)+!t  (i  +  »)+Mv;h. 


(B4) 


This  equation  also  has  a  root  u 


-  1  if  MV  *  MV‘  +  4  T  +  3T.  as  can  be 
o  Ao  J  eo  1 

easily  verified.  The  procedure  then  is  to  solve  Eq.  (B4)  numerically 

starting  from  this  root  and  following  it  as  VQ  increases.  Assuming  a 

singly  ionized  nitrogen  plasma  with  nQ  -  1015,  Tg0  =  2ev,  Tlo  =  0.5ev,  B  * 

2KG,  Fig.  B1  shows  a  plot  of  downstream  values  of  density  and  electron 

2  5 

temperature  versus  magnetosonic  Mach  number  m  ■  v/fv.  +  -=-  T  /M  + 


2 


3 


4tt 

Figure  B1  Downstream  electron  temperature  (dotted)  and  density  ratio 
(solid)  for  a  shock  with  =  5/3,  *  3  and  adiabatic  ion 

heating. 
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OICY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
OICY  ATTN  F.  REIBERT 


GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 
HMES 

COURT  STREET 
SYRACUSE,  NY  13201 
01CY  ATTN  G.  HILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AX  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
01CY  ATTN  J.M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFARD 
01CY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

01CY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DASIAC 
01CY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 
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MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CT  ATTN  N.  HARRIS 
Old  ATTN  J.  MOULE 
01CT  ATTN  GEORGE  MROZ 
Old  ATTN  W.  OLSON 
Old  ATTN  R.W.  HALPRIN 
Old  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
Old  ATTN  P.  FISCHER 
Old  ATTN  W.F.  CREVIER 
Old  ATTN  STEVEN  L.  CUTSCHE 
Old  ATTN  R.  BOGUSCH 
Old  ATTN  R.  HENDRICK 
Old  ATTN  RALPH  KILB 
Old  ATTN  DAVE  SOULE 
Old  ATTN  F.  FAJEN 
Old  ATTN  M.  SCHEIBE 
Old  ATTN  CONRAD  L.  LONGHIRE 
Old  ATTN  B.  WHITE 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NEW  MEXICO  87106 
01CY  R.  STELLINGWERF 
Old  'M.  ALME 
Old  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

Old  ATTN  JOHN  MORGANSTERN 
Old  ATTN  G.  HARDING 
Old  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  COR? 

12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01CY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01CY  ATTN  IRVING  L.  ROFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

Old  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
Old  ATTN  FORREST  GILMORE 
01CY  ATTN  WILLIAM  B.  WRICHT,  JR. 
01CY  ATTN  ROBERT  F.  LELEVIER 
01CY  ATTN  WILLIAM  J.  KARZAS 
01CY  ATTN  H.  ORY 
01CY  ATTN  C.  MACDONALD 
01CY  ATTN  R.  TURCO 
01CY  ATTN  L.  DeRAND 
01CY  ATTN  W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
Old  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01CY  ATTN  VINCE  TRAPANI 


57 


SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  PRIEMAN 
OICY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  COODRIDCE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 


OICY 

ATTN 

DONALD  NE1LS0N 

OICY 

ATTN 

ALAN  BURNS 

OICY 

ATTN 

G.  SMITH 

OICY 

ATTN 

R.  TSUNODA 

OICY 

ATTN 

DAVID  A.  JOHNSON 

OICY 

ATTN 

WALTER  G.  CHESNUT 

OICY 

ATTN 

CHARLES  L.  RINO 

OICY 

ATTN 

WALTER  JAYE 

OICY 

ATTN 

J.  VICKREY 

OICY 

ATTN 

RAY  L.  LEADABRAND 

OICY 

ATTN 

G.  CARPENTER 

OICY 

ATTN 

G.  PRICE 

OICY 

ATTN 

R.  LIVINGSTON 

OICY 

ATTN 

V.  GONZALES 

OICY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,  JR. 
OICY  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90273 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 
OICY  ATTN  D/  STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 


* 


